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ABSTRACT 
We examine the consistency of the thermodynamics of the most general class of con- 
formally flat solution with an irrotational perfect fluid source (the Stephani Universes). 
For the case when the isometry group has dimension r > 2, the Gibbs-Duhem relation 
is always integrable, but if r < 2 it is only integrable for the particular subclass (con- 
taining FRW cosmologies) characterized by r = 1 and by admitting a conformal motion 
parallel to the 4- velocity. We provide explicit forms of the state variables and equations 
of state linking them. These formal thermodynamic relations are determined up to an 
arbitrary function of time which reduces to the FRW scale factor in the FRW limit of the 
solutions. We show that a formal identification of this free parameter with a FRW scale 
factor determined by FRW dynamics leads to an unphysical temperature evolution law. 
If this parameter is not identified with a FRW scale factor, it is possible to find examples 
of solutions and formal equations of state complying with suitable energy conditions and 
O ■ reasonable asymptotic behavior and temperature laws. 

On 
O ! PACS numbers: 04.20.-q; 98.80.k 

cr- 

u 
bJQ. 

> 



On 



£ 



"Work supported by CONACYT, Mexico, project No. 3567-E 



I. Introduction 

The "Stephani Universes" is the generic name for a class of metrics comprising the 
most general conformally flat solution with an irrotational perfect fluid source 1-6 . These 
solutions, admitting (in general) no isometries, generalize FRW spacetimes (their particular 
case with vanishing 4-acceleration), and could have a physical interest as simple inhomoge- 
neous and anisotropic cosmological models. However, it is well known that Stephani Uni- 
verses (except for the FRW subcase) do not admit a barotropic equation of state p = p(p), 
where p and p are the pressure and matter-energy density, perhaps explaining why so few 
references are found in the literature 3 ' 7 ' 8 studying the physics of (non-FRW) Stephani Uni- 
verses. However, barotropic equations of state might be too restrictive 9 , and so we aim in 
this paper to verify whether one can find more arguments to gauge the physical viability 
of these solutions besides simply dismissing them for not admitting a barotropic equation 
of state. 

Coll and Ferrando 10 addressed the question of the consistency of the thermodynamical 
equations with Einstein field equations, deriving rigurously the criterion to verify if a 
single component perfect fluid source of a given exact solution admits what these authors 
denote a "thermodynamic scheme". However, Coll and Ferrando did not go beyond the 
admisibility of their consistency criterion, that is, into the physics of the fluids: note that 
it is perfectly possible to have unphysical fluids whose thermodynamics is formally correct. 
In a recent paper 11 , we have expanded and complemented the work of Coll and Ferrando 
by applying their criterion to irrotational, non-isentropic and geodesic fluids (the perfect 
fluid Szekeres solutions). Regarding the Stephani Universes, Bona and Coll 7 did apply 
the work of Coll and Ferrando to these solutions, claiming that non-barotropic cases only 
admit a thermodynamic scheme if r > 2 and presenting a very brief discussion of the 
thermodynamics of non-barotropic Stephani Universes with r = 2. By providing a specific 
counterexample, we prove in this paper that the result of Bona and Coll is incorrect (see 
section V) . We also provide a deeper discussion of the thermodynamics of non-barotropic 
Stephani Universes with r < 2 admitting a thermodynamic scheme. It is important to 
specify that the study of this thermodynamics assumes the matter source to be a single 



component perfect fluid. If the source were a mixture of perfect fluids, the study of 
its thermodynamics would involve looking at chemical potentials, and so would be an 
altogether different problem, a problem which will not be addressed in this paper. We 
assume henceforth that all mention of the term fluid indicates a single component fluid. 
The contents of this paper are described below. 

We present in section II a summary of the equations of the thermodynamics of a 
general relativistic perfect fluid, together with the conditions for admissibility of a ther- 
modynamic scheme. We re-phrase the conditions derived by Coll and Ferrando in terms of 
differential forms expanded in a coordinate basis adapted to the comoving frame in which 
the Stephani Universes are usualy described. These conditions are applied in section III 
to the Stephani Universes, yielding an interesting result: these solutions do not admit a 
thermodynamic scheme in general, that is, with unrestricted values of their free parame- 
ters. However, under suitable restrictions of these parameters, we find a specific subclass 
of non-barotropic Stephani Universes which does admit a thermodynamic scheme when 
r < 2. This subclass is the counter example to the work of Bona and Coll mentioned 
above. 

In section IV, we derive for the non-barotropic Stephani Univeres complying with a 
thermodynamic scheme explicit expressions of all state variables: p and p, particle number 
density n, specific entropy S and temperature T, as well as two-parameter equations of 
state linking them. The latter turn out to be difficult to interpret as there is no clue 
on how to fix the only time dependent free parameter of the solutions and its relation 
with the matter energy density. We explore in section V the strategy which consists in 
formally identifying these quantities with the scale factor and matter-energy density of a 
FRW limiting spacetime assumed to comply with a "gamma law" equation of state. The 
resulting temperature evolution law is unphysical and does not reduce in the FRW limit 
to that expected for a FRW cosmology with such an equation of state. On the other hand, 
if this identification with FRW parameters is abandoned, we show in section VI that it is 
possible to obtain temperature evolution laws and equations of state which, being more 
formal than physical, are not altogether unphysical and do comply with suitable energy 
conditions. Conclusions are presented and summarized in section VII. 



We prove in the Appendix that Stephani Universes admitting a thermodynamic scheme| 
are characterized invariantly by admitting: (a) a spacelike Killing vector (r = 1) and (b) a 
conformal Killing vector field parallel to the 4- velocity. In fact, we show that these Stephani 
Universes comprise the most general class of solutions with an irrotational perfect fluid 
source admitting this type of conformal symmetry. 

II. Thermodynamics of a non-isentropic irrotational perfect fluid. 

Consider the energy-momentum tensor for a perfect fluid 

T ab = (p + p)u a u b + pg ab (1) 

where p, p and u a are the matter-energy density, pressure and 4-velocity, respectively. 

ib 

;b 



This tensor satisfies the conservation law T a \ = which implies the contracted Bianchi 



identities 

p+(p + p)e = (2a) 

h b a p, b +(p + p)u a = Q (26) 

where © = u a a , u a = u a -^u b and h b a = 5 b + u a u b are respectively the expansion, 4- 
acceleration and projection tensor and p = u a p, a . The thermodynamics of a perfect fluid 
is essencially contained in the matter conservation law, the condition of vanishing entropy 
production and the Gibbs-Duhem relation. The first two are given by 

(nu a ). a = (3a) 

(nSu a ). a = (36) 

where n is the particle number density and S is the specific entropy. Condition (3a) 
inserted in (3b) leads to u a S, a = S = 0, so that S is conserved along the fluid lines but is 
not a universal constant. In the latter case we have: d5" = 0, and the fluid is isentropic, 



admitting a barotropic equation of state. The Gibbs-Duhem relation can be given as the 
1-form 



w = dS = f 



*(2)+"® 



(4) 



where T is the temperature. The necessary and sufficient condition for the integrability of 

(4) 

u> A do; = necessary and sufficient (5) 

subjected to fulfilment of the conservation laws (2) and (3), are the conditions which Coll 
and Ferrando denote admissibility of a "thermodynamic scheme" . These conditions were 
given by these authors as 

(pdp — pdp) A dp A dp = (6) 

Another integrability condition, not examined by Coll and Ferrando, is 

do; = sufficient (7) 

The perfect fluid source of the Stephani Universes is characterized by an irrotational 
(hence, hypersurface orthogonal) and shear-free 4-velocity. For such a fluid source there 
exist 2 ' 4 ' 5 local comoving coordinates (t, x l ), such that the metric, 4-velocity, 4-acceleration, 
expansion and projection tensor are given by 

ds 2 = -N 2 dt 2 + L 2 5 lJ dx l dx J (8a) 

0/3 

u a = N- 1 6? u a = (log N),Jt (86) 

hob = gijtfJl ( 8c ) 



where © = ©(£) and the metric function L is (in general) a function of all the coordinates 
(t, x % ). In this representation X = (l/N)X jt for all scalar functions and the Bianchi 
identities and conservation laws (2) and (3) become 



P,t + (p + P)Q-Ogn) 7 t = 



(9a) 



P, i + (p + p)(logN)^ = 



(96) 



n 



no(x l ) 
L 3 



(9c) 



S = S(x i 



(9d) 



where no(x l ) appearing in (9c) is an arbitrary function denoting the conserved particle 
number distribution. In the the coordinate basis of 1-forms (dt, dx l ) associated with the 
comoving frame (8), the Gibbs-Duhem relation reads 



u = S jdx 1 



T 



-) +»( i 



dx l 



(10) 



where the t component of to in this coordinate basis vanishes due to (9d). A sufficient 
integrability condition of (10) is given by 



du = W 1 



dt A dx i 



ti- 



nT 



+ Wi 



dx l A dx J 



«j 



nT 







(11) 



Wi 



ti 



P[,i n ,t] ~ n2T [,t s ,i] 



n 



/n A T t . n jt 

(P + P) — -T^r ~Ui — 
n 1 n 



T t , riA 
T n 



W, 



P[,i n ,j] + n2T [,i S ,J] _ T l,iP,j] 



ij 



~(P + P) 



T\, +Tuu\ n ,• 



n T \ T J n 

where square brackets denote antisymmetrization on the corresponding indices. The nec- 
essary and sufficient condition (5) is given by 
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^ dx l A dx j A dx k ^ dt A dx l A dx j 

dwAW = Xi» tffi + X ^ ^f2 = ( 12 ) 



X ti j = P[,tP,i n ,j] 

Conditions (12) are entirely equivalent to (6) provided by Coll and Ferrando. One can 
obtain the latter form the former simply by using (2) and (3) (in their forms (9)). However, 
(11) and (12) are more intuitive than (5) and (6), as they directly incorporate state variables 
such as n, S and T, and their relations with p and p. Condition (12) is also more practical 
than (6), as it is easier to use n and S from (9c) and (9d) than to compute the set 
(p, p, p, p) in exact solutions in which these quantities can be quite cumbersome. The 
sufficient condition (11), not examined by Coll and Ferrando, is also helpful, since if its 
fulfilment guarantees that (6) (or (12)) holds. 

As shown in the following section, if a solution of Einstein equations is available (thus 
providing p and p in terms of the metric functions) it is straightforward to verify the 
admissibility of the thermodynamic scheme. This we will do for the Stephani Universes 
which are particular cases of (8), and to do so we suggest the following procedure: (a) 
solve the conditions (12) and substitute the solution into (10), thus identifying possible 
(non-unique) forms for S and T; insert the obtained forms of T and n into (11) in order 
to verify if further restrictions follow from the sufficient conditions. If these conditions 
hold, the equations of state linking the state variables (p, p, n, S, T) (together with their 
functional relation with respect to the metric functions) follow directly from integrating 
them. 

III. The Stephani Universes. 

The Stephani Universes are described by the particular case of (8a) given by: 

L(t ' X ^ = 1 + 2A i x i + (A 2 + (k/A)E 2 )5 lJ x^x^ (13) 
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with Ai(t) = (A x (t),A y (t),A z (t)), A 2 = 5 ij AiAj, k(t) and R(t) are arbitrary functions. 
Notice that the Stephani Universes contain FRW spacetimes as the particular case A, = 0, 
k = ko = const, in (8) and (13). As it is well known that the latter are the subclass of 
barotropic Stephani Universes, we assume hereafter (and unless stated otherwise) that all 
mention of Stephani Universes excludes their FRW subclass. The field equations associated 
with (8) and (13) are 



e 2 

p(t) = ^ + 3k 



(14a) 



P 



~P 



P,t 



(146) 



3L, t /L 
where the contracted Bianchi identity (9a) has been used as a definition of p. The remaining 

Bianchi identities and conservation laws are given by equations (9b-d), while the Gibbs- 

Duhem 1-form (10) and the integrability conditions (11) and (12) follow as 



sv 



(p+p) fl 



n 



(15) 



W t/ , 



(p + p)^- -p, t 



n,i-{p + p)Uin, t = 



(16a) 



where 



Wij = 0^O[ l (logn) )j] =0 



(166) 



^ = (logT) ;i + ^ = (logTAr) ; 



(16c) 



X t 







. tij = v =? p[,in tj ] = => (logiV)[ ;i (logn)j] = (17) 

Inserting the forms of iV and n given by (8a) (9c) and (13) into the necessary and sufficient 
condition (17) yields a general solution of the latter given by 



a(t) lo£ 



tI+wMi? 



log(c(t)) 



(18) 



where (a, 6, c) are arbitrary functions. This integrability condition is not satisfied in gen- 
eral, that is for arbitrary forms of the free functions (Ai,k,R) appearing in the metric 
function (13). Bona and Coll 7 have claimed that this condition is only satisfied by Stephani 
Universes with isometry groups of dimension r > 2. However, it is possible to provide a 
particular case of this metric (characterized by r = 1, see Appendix) which satisfies (16) 
and (17) and so, leads to well defined forms for T and S. This case is characterized by the 
existence of a conformal Killing vector field parallel to the 4- velocity (see Appendix), the 
corresponding forms of Ai and k are 



ko_ bo 
R 2 R 



A % = a z R, k = -£. + JL - 4d ij a,aj (19a) 



where (a«, ko, 6o) are arbitrary constants and R remains arbitrary. With these parameter 
values, (18) holds with 



a 



(t) = 1, b{t) = 1/3, c(t) = ~{1/R),t (196) 



"-o = / 3 where: / = 1 + \k 5ijX l x j (19c) 

Inserting (19a) and (19c) into (15) yields 

T^ = -3(p + p)(L/) 4 (F//) ;i (20a) 

where 



F = 2a i x i + lb 8 ij x i x j (206) 

This equation shows how non unique forms of S and T emerge if we demand only the ful- 
filment of condition (17) (or (18)). A general expression for these quantities, in agreement 
with (15), is given by 

S = S(a) a = -j (21a) 

9 



where a prime denotes derivative with respect to a and (9c) (13) have been used to eliminate 
L in terms of n. Regarding the sufficient conditions (16), T given by (21b) satisfies (16a). 
This is easily verified by eliminating iii from (16c) and inserting this result together with 
Tj/T computed from (21b) into (16a). On the other hand, inserting (21b) into (16b) and 
using (17) leads to the condition S" = 0, so that S is a linear function of a. This yields 
the following forms for T and S 

S = S + a (23a) 

T = «E+p. (23.) 

which are compatible with both sets of conditions (16) and (17). A discussion on the inter- 
pretation of the expressions derived above is provided in the following sections. 

IV. Formal equations of state. 

The particular subclass of Stephani Universes complying with the thermodynamic 
scheme, as characterized by the parameter restrictions (19a), is described by the confor- 
mally FRW metric 



ds 2 = ^ 



-dt 2 + 



(1+JWy^Xj) 2 



(24a) 



1 + jkoSijX l x J f 



$ = , 4 " ; r = (246) 

l + \koSi j x % x3 + R(2a i x % + \boSi j x % x3) f + RF 

Its corresponding field equations are 



■?. 



b>=(it) +| + |-4^a, (25a) 



p = -p-lRp, R (l-Ra) (256) 
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where we have chosen 0/3 = Rt/R and a is given by (21a). These Stephani Universes can 
be characterized invariantly by the existence of a conformal Killing vector field parallel 
to the 4- velocity (see Appendix). Assuming the forms of T and S from (23), the term a 
appearing in (25b) (and in other expressions) is a linear function of S. Also, since p and 
R are both functions of t, all terms involving R and p^n can be expressed as functions of 
p. This results in the following forms for generic equations of state, expressing p, n and T 
in terms of p and S 

p(p, S) = -p- \Rp, R [1 + (S - S )R] (26a) 



n(p, S) 



i + S-So' 



(266) 



T(p, S) = ~ 3R5p > R 3 (26c) 

These are formal equations of state, in the sense that they are not dictated by physical 
considerations and imposed before solving the constraints of the field equations (physical 
equations of state), but arise as a consequence of imposing the fulfilment of the thermody- 
namic scheme on metric functions whose spacial dependence has been fixed by imposing 
conformal flatness (a geometric constraint: vanishing of the Weyl tensor). The best one 
can do in this case is to verify if these formal thermodynamic relations could be manipu- 
lated in such a way that solutions (24) could describe physically reasonable cosmologies. 
However, equations (25) are still undetermined: a choice of p(R) must be made in order to 
determine these equations and to be able to integrate the Friedmann-like equation (25a) 
to yield the time evolution of the metric. Unfortunately, there is no clear cut way guiding 
one on how to select p(R). Various possibilities are explored below. 

V. FRW limit. 

The metric (24) bears a close resemblance to a FRW metric. As equations (24)- 
(25) reveal, the term a, related to the spacially dependent entropy density, is the term 
that makes these solutions inhomogeneous and anisotropic (i.e. non-FRW). In fact, their 
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FRW limit follows if 60 — ^ and ai —> 0, so that a — > 0. Under this limit, the metric 
and field equations (25) become the metric and field equations of a perfect fluid FRW 
spacetime with scale factor R and ko = 0, ±1 marking the curvature of the spacial sections. 
This correspondence and resemblance to FRW cosmologies motivates us to verify if these 
solutions could be considered as some sort of "near-FRW" cosmologies and if they could be 
examined within the framework of a FRW limit. Hence we pose the question of whether 
the state variables and their equations of state become those one would expect of a "near- 
FRW" cosmology if we assume that R can be fixed as it were a FRW scale factor. Along 
these lines, we notice that the state variables p, p and n, given by (25a), (25b), (26a) and 
(26b), also tend to their FRW values, but the limiting form of T in (26c) takes the strange 
form T — > — R 5 p #. In order to see what sort of temperature law and equations of state 
correspond to this FRW limit, we assume a "gamma law" equation of state po = (7 — l)p, 
where po is given by setting S = Sq in (26a), the form of p(R) becomes 

R^~ 31 



and so the various forms (26) of the equation of state become 

p(p, S) = (7 - l)p + § 7 #op 1 - 1/37 (S - So) (27a) 



n(p, S) = Ip 1 ^ + S-So} 3 (276) 



d4 I-I/37 

T(p, S) = —^ 3 (27c) 

Irrespective of the interpretation of these strange formal thermodynamic relations, notice 
that p and n do reduce to their FRW values in the FRW limit, though T — >■ p 1-4 / 3 ')' oc 
.R 4-37 , a temperature evolution law which has no relation to that expected at the FRW 
limit: for 7 = 1 (dust), T oc R instead of T = const, and for 7 = 4/3 (radiation), 
T oc const, instead of T oc R~ l . Therefore, this formal identification of R with the FRW 
scale factor leads to contradictory results. Of course, one could also consider R as a FRW 
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scale factor associated with other equations of state, however, the presence of the strange 
term —R 5 pji in (26c) makes it highly improbable for this temperature evolution law to 
have any meaningful correspondence with that of a limiting FRW cosmology. 

VI. Examples complying with minimal physical requirements. 

Looking at the equations of state (26) as formal thermodynamic relations, which 
hopefuly might provide at least a gross approximation to physical relations, and if one is 
prepared to abandon the identification of R with a specific FRW scale factor, it is possible 
at least to verify if the free parameters (a*, ko, bo) and R(t) can be selected to assure that 
solutions (24) and equations (26) comply with minimal physical conditions, such as energy 
conditions and an acceptable asymptotic behavior. In order to examine equations (24)- (26) 
within this framework, consider the coordinate transformation 

x = 2r(x/2) sin 9 cos <p 

s = 2r( x /2)si„*si„„ r(x /2) = {££$, *£=-! 

z = 2r( X /2)cos6 
which brings the metric (24) into the simple form 



-dt 2 + R 2 [d X 2 + T?{x){d9 2 + sin 2 9d<p 



oV = L — -^ -^- (28a) 

[l + R(2E(x)V(e,^ + b ^(xm)} 



E w=te k i:=-i (286) 



V(6, ip) = a\ sin 9 cos cp + Cb2 sin 9 sin <p + as cos 9 (28c) 

The state variables in (26) are unaffected by this coordinate transformation, with a 
and the function n$ in (19c) which provides the conserved particle number at an inicial 
hyper surface t = const, now given by 

a = 2E( X )F(0, ip) + 6 S 2 (x/2) (29a) 
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"° = {co S h«(I X ) *, = -! < 2% » 

From (28c), the term V containing the dependence on (9, <p) is bounded, hence a = —F/f 
might only diverge along x ~^ =t°° (irrespective of the value of p — pit)) for the case 
ko = — 1. This fact can be examined from another angle: the function a (and so, the 
entropy density S) can be expressed in terms of no as 

S = S + a = S + 4nJ /6 (l - /e ™J /3 ) 1/2 + b (l - /e nj /3 ) (29c) 

and so, for ko = —1, the initial particle number no becomes infinite as x -* ±oo, this 
infinite concentration of particles causes the entropy density S and the pressure to diverge 
along these limits. This means that solutions with ko = —1 have an undesirable asymptotic 
behavior. On the other hand, the locus R = marks another singularity, analogous to a 
FRW big bang, while the vanishing of the denominator in (28a) indicates an asymptotically 
deSitter evolution (p — ► — p, n — ► oo) characterized by the unphysical behavior T —> oo. 
Therefore, for the formal equations of state (26) to have any physical meaning, we must 
choose ko — 1 and demand the condition 1 — Ra ^ to hold, together with the dominant 
and weak energy conditions which can be combined into the restriction: < p/p < 1. 
Also, all state variables must diverge at the big-bang singularity R = 0. From equations 
(26), the conditions < p/p < 1 and T —> oo hold at the limit R — * if p ~ i^ _ ( 4 + m ) for 
?n > at this limit, leading to the following asymptotic values: 



-»|(m+l) (30a) 



n » _L « p 3/(m+4) (306) 



T^^i^(m + 4)^/(-+ 4 ) (30c) 

as R — > 0. Since R is no longer constrained to be interpreted as a sort of FRW scale factor 
in a FRW limit, we can devise a simple example complying with the conditions mentioned 
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above and the asymptotic limits (30) by choosing a simple power law p = (Rq/R) 6 , where 
Ro is a constant. This choice leads to the following state variables 

p(p,S)=p + 2R Q (S-So)p 5/6 (31a) 



n(p, S) = [p 1 / 6 + S - So] 3 (316) 



T(p,S)= r - °: -3 (31c) 



6^ p 5/6 
[pl/6 + 5 _ 5o ]« 

which yield a stiff fluid equation of state p/p — > 1 and T — > oo near the big-bang singularity 
(i? = and/or p — > oo). However, as R —> oo (or p — ► 0), the dominant energy condition 
could be violated: p/p — > oo. This can be avoided by selecting the remaining arbitrary 
constants (a^&o) m such a way that So > S and 1 — i?a = 1 + it^o 1 — So) > holds 
everywhere and the Friedmann-like equation 

(2f) 2 = 4L - 4 - I + 4^^ (32a) 

\ R J 3R G R 2 R J y ' 

obtained by substituting ko = 1 and p = (Rq/R) 6 into (25a), has no solutions _R(£) allowing 

for i? — > oo, or equivalently, that p does not vanish along the time evolution of the fluid. 

These conditions require (32a) to have a real positive root and 

% = -" (^-X - "^ + ^Wi < (326) 

R 3 \ R J 2R J y ' 

along this root, hence R(t) is convex. If we assume R = Rq > to be the value along 

which Rj/R = 0/3 vanishes, equation (32a) fixes Rq in terms of the parameters (aj,6o) 

as the positive root of: 

(l + 45 lJ a l a j )R 2 l -b R -l = (32c) 

which substituted into (32b) yields the condition of convexity. It is not difficult to find 
combinations of parameters (a^,6o) so that the fluid has the desired type of kinematic 
evolution (qualitatively analogous to a standard "closed" FRW cosmology) and physically 
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correct behavior of the state variables: that is, to have the ratio p/p = 1 and T — ► oo at 
the big bang evolving to values in the range < p/p < 1 and cooling to T finite as the fluid 
reaches its maximum expansion, bounces and then recollapses with p/p = 1 and T — »■ oo 
at the big crunch. Such an example is illustrated in figure 1. 

VII. Conclusions. 

We have investigated the consistency of the thermodynamic equations following from 
the condition of existence of a thermodynamic scheme for the Stephani Universes, whose 
source is a non-isentropic perfect fluid (thus, not admitting a barotropic equation of state). 
This work has aimed at improving the study of this type of solutions, as classical fluid 
models generalizing FRW cosmologies, in contrast to a widespread attitude of simply 
disregarding them for not admitting a barotropic equation of state. 

For the particular subclass of Stephani Universes admitting a thermodynamic scheme, 
the resulting equations of state have an ellusive interpretation, as there is no blue print 
on how to select the time dependent free parameter of the solutions, an arbitrary function 
R reducing to the FRW scale factor in the FRW limit. We have shown that by formally 
identifying this parameter with the FRW scale factor of a FRW cosmology satisfying a 
"gamma law" equation of state leads to an unphysical temperature evolution law, totally 
unrelated to that of their limiting FRW cosmology. The question of how to select these 
parameters in a convenient way remains unsolved, though the adequate theoretical frame- 
work to carry this task has been presented in section VI. We have shown that combinations 
of free parameters exist so that the formal equations of state comply with minimal physical 
requirements. 

We have shown that Stephani Universes (other than the FRW subclass or the sub- 
class presented in previous sections) are not compatible, in general, with a thermodynamic 
scheme. This fact seems to disqualify these solutions as classical fluids of physical interest. 
However, the latter can still be useful if they are examined under a less restrictive frame- 
work than that of the simple perfect fluid. In this context, the Stephani Universes (like 
other perfect fluid solutions with a shear-free 4-velocity) can be recast as exact solutions 
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for a fluid with a bulk viscous stress 12 . From this point of view, the thermodynamics is not 
only totally different as that of the perfect fluid case but much less restrictive and more 
amenable to satisfy the criteria for constructing inhomogeneous and anisotropic cosmolog- 
ical models of physical interest. 
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Appendix. Isometries and conformal symmetries, 
(a) The metric (24) admits a one-parameter group of isometries. 

Being conformally flat, (24) admits a Gi$ of conformal symmetries. Various isometry 
subgroups are readily identified. If ko = 0, irrespective of the values of the remaining 
constants a,i, (24) becomes spherically symmetric (this can be verified through re-scalings 
of the form x % = x l + c l ). Consider now (24) with ko = ±1 and Oj arbitrary nonzero 
constants. Since we are assuming this spacetime to be non-static (© = 3R jt /R ^ 0), its 
Killing vectors (if they exist) must all be tangent to the hypersurfaces orthogonal to the 
4-velocity 15 . Spherical symmetry is easily identified by setting a^ = and bo 7^ 0, while 
FRW subcases follow from a^ = = bo. However, for general values of these constant 
parameters, it is not easy to identify at first glance the existence of isometries groups G\ 
or G3 acting along the hypersurface orthogonal to the 4-velocity. 

In order to find out if (24) admits isometries, and if so, the dimension of their orbits, 
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we follow the result of theorem 3 in Bona and Coll 6 . Identifying carefully the parameters 
used by these authors: a, b and 4> (see their equations (2)-(4)) with the corresponding 
parameters in (8) and (13), we find that their function a is our function R, their vector 
b is the vector formed by the functions Ai in (13), while = A 2 + (k/4) * R 2 , that is, 
the coefficient of the quadratic term in the denominator of (13). The case complying with 
the thermodynamic scheme (the metric (24)), is defined by (19a), thus we have for this 
metric b = (ai,a2,a,3)R and (f> = boR. Therefore, metric (24) with arbitrary values of its 
free parameters R, a^, bo, corresponds to the case: b ^ 0, b ^ with b A b = and with 
4> = boR y^ 0. According to the classification of page 616 of Bona and Coll 6 , the metrics 
(24) admit a one parameter isometry group with r = 1, a case these authors identify as 
axially symmetric. 

(b) All Stephani Universes associated with the metric (24) admit a conformal 
Killing vector parallel to the 4-velocity. 

We have proven that the metric (24) describes a class of Stephani Universes admitting 
a thermodynamic scheme and having an isometry group of dimension r < 2. On the other 
hand, it is well known 16 that the 4-velocity in FRW spacetimes (u a = Sf in the comoving 
coordinates of (24)) is a conformal Killing vector, satisfying £(a;6) = V\ )#a6 with scale 
factor ip (0) = R t . Since (24) is conformally related to a FRW metric, the vector field 
£ a = 8f in (24) (parallel to the 4-velocity u a = (— <&)~ l 5f) is a conformal Killing vector £ a 
in (24), satisfying £( a; &) = ipQab, with conformal factor i/j given by 16 



i> = v%) +e a (iog*), a = Rt (i - j^rf) ( A1 ) 



(c) The metric (24) describes the most general perfect fluid spacetime admit- 
ting a conformal Killing vector parallel to the 4-velocity. 

It is known 17 that the existence of a conformal symmetry of this type ( £ a = 0,u a , 
satisfying £( a; b) = ipQab) requires the fluid 4-velocity to be shear-free, with the remaining 
kinematic parameters given by 
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u a = u a ., b u h = h h a {\ogQ), h (A2) 

ab = U[a;6] + U[ a U h ] = 0(0" 1 ^[ a ) ;6 ] (A3) 

ie = i^ a = ^(logO), a = ^ (A4) 

Consider the irrotational case u a b = 0, as described by the comoving coordinates of (8) 
but with L not necessarily equal to (13). The specific form of the 4-acceleration in (A2) 
implies the constraint 

(i/L) it = jeflj(x') 

where J is an arbitrary function. Inserting this constraint into (8), and demanding the 
fulfilment of the Einstein field equations Gij = and Gu — Gjj = (i ^ j), the metric (8) 
becomes (24a) with J = f. 
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Figure Caption 

Figure 1. Thermodynamic quantities associated with the example presented 
in section VI. 

We chose the following values of the constant free parameters: a\ = 10 -4 , a^ = 
9 X 10" 5 , a 3 = 1.1 X 10" 4 , b = -9/20, so that the positive solution of (32c) is #0 « 1.18. 
Figure (la) displays 0/3 = Rj/R vs. R, showing that the latter is bounded between 
R — and R = Rq, with the latter constant being the value of R at which © vanishes. 
This means that p evolves between p — 1 at R = Rq and infinity at R = (big bang 
and big crunch singularities). Figure (lb) exhibits the ratio p/p in terms of R and the 
"radial" coordinate x, for 9 = 7r/2 and = 0. Since the constants o» are very small, the 
plot looks qualitatively analogous for all values of these "angular" coordinates. Notice how 
the constraint < p/p < 1 holds throughout the evolution of the fluid, from p/p = 1 
at R = 0, keeping the same value along x = and evolving to p/p ~ at R = Rq and 
X ~ ±7r. Figures (lc) and (Id) display equations of state (31a) and (31c) The ratio p/p 
and T are plotted in terms of arctan(p) in the range (arctan(l), ir/2) for various values 
of S = Sq — a in the range between S = and S — So — 9/20. Notice that p/p — ► 1 
and T — > 00 as arctan(p) — »■ 7r/2 (that is: p —> 00) for all values of S. As p — »■ 1, T and 
p/p decrease at various rates depending on the value of S. The parameter values in these 
plots do not follow from any physical consideration, we simply aim to illustrate that free 
parameters exist so that the solutions satisfy the minimal physical conditions discussed in 
section VI. These plots were obtained with the symbolic computing program MAPLE. 
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